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ABSTRACT 



We consider the U ( 1 ) charged sector of ABJM theory at finite temperature, which corresponds 
to the Reissner-Nordstrom AdS black hole in the dual type IIA supergravity description. In- 
cluding back-reaction to the bulk geometry, we show that phase transitions occur to a broken 
phase where SU(4) R-symmetry of the field theory is broken spontaneously by the condensa- 
tion of dimension one or two operators. We construct the composite operators out of fields in 
ABJM theory and describe the phase transition with the dual gravity solutions. We show nu- 
merically and analytically that the relevant critical exponents for the dimension one operator 
agree precisely with those of mean field theory in the strongly coupled regime of the large N 
planar limit. 



1 Introduction 



The ABJM theory is the three dimensional 5\£ = 6 U (N) x U(N) superconformal Chern-Simons 
theory with level (k, —k) and dual to the type IIA string theory on AdS4 x CP3 background HI. 
Some test of this duality has been carried out largely based on the integrability with indication 
of some additional structures compared to the well known AdSs/CF^ counterpart O-im. The 
type IIA supergravity description is dual to the large N planar limit of the Chern-Simons theory 
where one takes N, k — > °° while holding 't Hooft coupling X = N/k fixed. 

Some probes of the Chern-Simons plasma at finite temperature are carried out recently via 
the consistent, CP3 invariant dimensional reduction of the type IIA supergravity 10. Alter- 
natively a finite temperature plasma can be completely characterized by towers of static length 
(1/mass) scales. They are arising as decaying spatial length scales of a perturbation theory when 
local operators are inserted at a certain point of the plasma. In Ref. [7], these scales are scanned 
for the low-lying bosonic modes from which the true mass gap m g (the lowest in all) and the 
Debye screening mass (the lowest in CT odd sector of the theory j8]|) are found for the 
Chern-Simons plasma. Including Yang-Mills plasmas, these two scales are well representing 
the universal characteristics of a certain strongly coupled plasma. For instance, the ratio mo/nig 
for the $1 = 4 SYM theory and the two-flavor model of QCD are matching with each other in 
the strong coupling limit supporting such a picture [|9]|. Leading thermodynamic corrections of 
ABJM theory in the small X expansion are explored in Ref. ifTOll . 

In this paper, we note that the ABJM theory possesses SU(4) xU(l) R-symmetries and 
consider its particular sector in which one turns on a finite U ( 1 ) number density. For the type IIA 
supergravity side, this sector is described by the U(l) charged Reissner-Nordstrom (RN) AdS 
black brane solution whose physics will be our basic concern for the study of the field theory at 
a strong coupling [6] . We begin with the thermodynamic stability of the RN AdS black brane 
solution and prove that it is thermodynamically stable at all temperatures including the zero 
temperature limit. This is contrasted with the R-charged black holes in the type IIB supergravity 
which is dual to the R-charged sector of the 9*1 = 4 super- Yang-Mills (SYM) theory. The R- 
charged solution becomes thermodynamically unstable below a certain temperature and the 
validity of the solution itself will be lost completely [fTTT|. Thus the RN AdS solution of this 
note is, up to now, the only known example where one has the thermodynamic stability of the 
gravity solution at all temperatures and the dual field theory description is precisely known at 
the same time. In other words, we found that the RN AdS black brane solution is relevant 
dual gravity description for the finite temperature ABJM theory with a chemical potential in the 
whole temperature region including the zero temperature, which is not the case for the 5\£ = 4 
SYM theory. 

We go on to study phase structure and transitions occurring within this charged sector. These 
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phase transitions turn out to be a second order type where the SU (4) R-symmetry is broken 
spontaneously lfT2ll . For this we scan behaviors of supergravity modes which are dual to the 
primary operators of some dimension A in the RN black hole background ifTSI . Among them, 
one finds the bosonic modes with mass squared ranged over — 9/4 < m 2 < — 3/2 which are 
responsible for the phase transition. As will be further explained later on, the only possibility 
for the present case is m 2 = —2 corresponding to one A = 1 operator of SU (4) representation 15 
and two A = 2 operators of SU (4) representations 15 and 84, whose detailed operator contents 
will be specified below. The transition related to the condensate of the A = 1 operator occurs at 
a critical temperature higher than that of the A = 2 operators. Above the transition temperature, 
the scalar field has to be set to zero to satisfy the required boundary condition while the other 
part of the original RN black brane solution remains intact. Below the transition temperature 
the scalar field begins to develop a profile whose boundary behavior represents a condensation 
of operator expectation value without introducing any external source field. Of course the ge- 
ometry will then be back-reacted accordingly. The solution represents the condensation of the 
operator expectation values below the transition temperature, which takes a particular direction 
in the space of the SU (4) representation 15. Thus the SU (4) symmetry is spontaneously broken 
down to its little group. We then study the critical exponents of the phase transitions and show 
both analytically and numerically that the exponents precisely match with those of the mean 
field theory. We shall also show that there is further symmetry breaking phase transition of the 
same nature at a lower temperature which involves condensation of the A = 2 operators. 

In Section 2, we shall discuss some relevant properties of the ABJM theory. Section 3 deals 
with the gravity description of the RN black brane and related physics in the field theory. In 
Section 4, our system described by RN black brane solution is shown to be thermodynami- 
cally stable. In Section 5, we discuss the development of gravitational instabilities of scalar 
modes below the critical temperature. This will fix the critical temperature which depends on 
the dimension of the corresponding operators. In Section 6, we discuss the phase transition 
by studying the gravity solution representing the condensation of operator expectation values 
including the back reaction to the bulk geometry. The critical exponents are shown to agree 
precisely with those of the mean field theory. Last section is devoted to the interpretations and 
concluding remarks. 

2 Thermodynamics of ABJM theory at small X 

The on-shell degrees of the ABJM theory consist of bosonic and fermionic matter fields Y 1 
and *¥] (I = 1,2,3,4) together with two gauge fields A m and A m . The complex scalar fields 
Y 1 are in the representation of (N, N, 4) of the U(N) xU (N) gauge as well as the SU (4) R- 
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symmetries. There are also the complexified Majorana fermions which are in the represen- 
tation of (N, N, 4) . The gauge fields A m and A m in the adjoint representations of the first U (N) 
and the second U (N) respectively, are coupled to the matter fields <I> = (Y 1 , ¥/) by 

D m f = d m <$> + iA m * - i®A m . (2.1) 

For further details of the Lagrangian and notation used in this note, see for example Ref. |f3l. 
The theory possesses global 9\£ = 6 3d superconfomal symmetries of OSp(6\4) whose bosonic 
part is given by the 3d conformal symmetry of 50(3, 2) multiplied by SU (4) R-symmetry. This 
corresponds to the isometry of the AdS4 x CP3 in the type IIA supergravity side. But one 
crucial difference from the 5\£ = 4 case is the fact that there is an extra global U(l) symmetry. 
The relevant charge is associated with the U(\) phase transformation of complex fields by a 
overall phase factor. Denoting overall U(l) parts of gauge fields by A^^ and A^ 1 ) the above 
charge is only coupled to the relative U(l) gauge field 

Am A m . 

Now note that the fields Y 1 and *F/ carry the U(l) charges (— 1,+1) or —1 in terms of 
the relative one. Due to the Gauss law constraint of the Chern-Simons theory, these charges 



are always accompanied by U(l) magnetic fluxes^! (+1, — 1). Thus for general k including 
even the nonabelian contributions, the basic degrees in the field theory in the deconfined phase 
exhibit an anyonic nature due to the (generically nonabelian) statistical interactions between 
them. But due to the large N planar limit where we send iV and k to infinity at the same time, the 
statistical interactions drop out since they are of higher order in l/k. Therefore, for instance, 
the total effective number of degrees in the weakly coupled small X limit is simply proportional 
to N 2 . Namely the entropy density at temperature T takes a value [JTJ 

s (1^0) = 2 ^1n 2 T 2 (2.2) 
71 

in the X — > limit. (The free energy density is always related to the entropy density by J = 
—S T/3 as dictated by the conformal symmetry.) On the other hand in the strongly coupled 
large X region, the system is described by the black brane solution in the gravity side whose 
Bekenstein-Hawking entropy density reads [OQ 

5W = 2^^1. (2 .3) 

v ; 21V2 VI 

The appearance of the suppression factor is not understood from the direct computation of 
the field theory. This drastic change of the number of degrees might be related to some remnant 
of anyonic interactions but we do not have any supporting evidence for this picture. 

Our main concern in this paper is the sector of the ABJM theory where one turns on a finite 
U(l) number density or equivalently the corresponding chemical potential /j. First consider the 



'The unit flux for the current case with e = 1 is given by h/e = 2nh = 1% where we set H = 1. 
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small X region of planar limit and the high enough temperature with T ^> p. The fermions can 
be ignored in this limit. In the thermal circle compactified effective theory, the fermions can 
be integrated out at weak coupling and high temperature limit as their Matsubara frequency 
starts with nT. Of course the U(l) current gets contributions from both bosonic and fermionic 
degrees of the theory. Their currents are conserved not separately but only in sum. Hence at 
weak coupling of small A,, instead of building up a fermi surface, occupation of bosonic states 
is energetically preferred. The scalar fields in this effective 2d description acquire a mass |fl4l . 

m 2 (r) = V + '«r. (2-4) 
where rrij is the thermal mass correction. For X <C 1, the thermal mass has the expression [[Toll 

m 2 = 11 8 _ (log X) 2 T 2 + 0(X 2 log X). (2.5) 

The theory lies in the unbroken phase if m 2 (T) > 0. On the other hand the system in the 
symmetric phase becomes unstable when m 2 (T) < Oor/j <C T < \J sjg ^ ■ As argued in 
Ref. Ifl4l . some of the operator of the field theory may acquire nonvanishing expectation values 
possibly leading to SU (4) R-symmetry broken phase. But the precise fate of the system at weak 
coupling requires a further study, which is beyond scope of the present paper. 



3 RN black brane and type IIA supergravity on AdS4 x CP3 

In the strongly coupled region of large X, the description of the ABJM theory by the type IIA 
supergravity on AdS4 x CP3 is appropriate since geometry there is weakly curved. The type 
IIA spetra compactified on CP3 space was known quite some time ago [fT3ll . Each mode of the 
resulting 4d supergravity is dual to a gauge invariant primary operator whose scaling dimension 
A is protected against quantum corrections. The presence of these bulk modes shows how the 
basic degrees of freedom are organized in the strongly coupled side of the ABJM theory. Thus 
study of these bulk modes for a given supergravity background will be our main tool probing 
the physics at the strong coupling. 

The bulk modes consist of infinite towers of spectra from spin zero to spin two. We shall be 
here briefly describing some of relevant low-lying modes for later discussions. Let us begin with 
the case of spin one: The lowest are two massless bulk gauge fields that are dual to the A = 2 
current operators in the field theory side. One is for the current of SU (4) singlet representation 
[ (000) in the SU(4) Dynkin label notation], which is identified with that of the extra U(l) 
global symmetry. In the gravity side the U(l) is related to the U(l) isometry of the M-theory 
circle from the lid perspective and its charge is carried by DO branes in a rough sense. As was 
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shown explicitly in Ref. [6], this gauge field is arising as a linear combination of 



■A™ + 3A°\. 



(3.1) 



where and A® 4 respectively coupled to the DO branes and DA branes wrapping CP2 four 
cycle inside CP3. The other combination A = A® —A® 4 becomes massive by Higgs mechanism 
with m 2 = 12, which couples to the A = 5 boundary current operator. It should be also noted 
that the monopole operator with overall field theory U(l) charges n(k, —k) (n E Z) is the one 
example of heavy BPS state coupled to this U(l) bulk gauge field JTJ. The second massless 
gauge field is in the adjoint 15 [(101)] of SU (4) and couples to the boundary SU(4) current 
operator of dimension A = 2. 
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Table 1: The low lying spectra up to the operator dimension 3 are presented. The upper and lower 
indices denote respectively the parity and the mass squared value m 2 of the supergravity mode. The 
whole spectra of CP3 singlet sector are presented in addition. 

For the bulk scalar modes, the lowest ones with m 2 = —2 are relevant for our later discus- 
sions, which correspond to A = 1, 2 operators. The A = 1 scalar mode couples to the boundary 
operator of SU (4) representation 15 [(101)] which takes the form 



tfj = Tr Y I Y] - (trace part) 



(3.2) 
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and has 15 independent components. The A = 2 modes involve two field theory operators of 
SU (4) representation 15 [(101)] and 84 [(202)] whose operator contents read 

0\ = Tr - (trace part) , 

d[ L = Tr Y {I Y^Y J ^ YX - (trace part) (3.3) 

where the trace part denotes any contractions between the upper and the lower indices. 

The fermionic bulk modes start with \m\ = corresponding to the operator dimension A = §: 
These are boundary operators of SU (4) representations 10 [(002)], 6 [(010)], and 10 [(200)]. 
The low lying modes up to A = 3 are listed in Tabled] 

A few comments are in order. We note that there are no bulk supergravity modes that are 
charged under the massless U ( 1 ) gauge field. The knowledge about the compactification spectra 
tells us about only linearized fluctuations of modes above the AdS4 or the AdS4 black brane 
solution. In Ref. [6], a consistent CP3 compactification keeping all SU(4) invariant modes is 
carried out explicitly. Any solutions of this 4d system can be consistently embedded into the 
original lOd supergravity theory. 

Our starting Lagrangian for the further discussion is 

L = 2F (* + 6 " - E ( < V *«) 2 + »<5 ^ ) ) (3-4) 



( /=l 



1T*= * . (3.5) 



when 



6tiV2X 

In this action, the Einstein Maxwell part with the negative cosmological constant is a fully 
consistent truncation [6] while the remaining scalar part is only valid up to quadratic order. 
Below we shall show that, above the critical temperature of the A = 1 operator, the stable 
solution of the system is given by the RN black brane solution with vanishing scalar fields. 
This part of the solution is fully consistent as just stated. Below the critical temperature, the 
relevant scalar field begins to develop and the corresponding set of solutions is valid only if the 
magnitude of the scalar field is small enough. But the set of solutions in the near critical region 
carries all the information about universal natures of the phase transition including relevant 
critical exponents. 

Similar models have been discussed many times in the bottom up approach lfT5l [T6l l22l . 
However, it is hard to identify the dual field theory and the corresponding operator for conden- 
sation in this bottom up approach. In our work, on the contrary, the identification of the physics 
of boundary CFT is straightforward. 



2 In the bulk, any dimensionful quantity is measured with respect to the AdS radius scale I which we set to be 
unity for the notational simplicity. 
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Below we take the following ansatz, 

ds 2 = e 2A ^(-h(r)dt 2 + dx 2 + dy 2 ) + ^-, 
\ ) h(r) 

A t = A t (r) , 4> a = 4> a (r), (3.6) 

to describe a finite temperature system with plane plus time (M 2 x W) translational symmetries. 
Plugging the above ansatz into the equations of motion, we are led to the following equations 
E3: 

z i=i 

h" + 3A'h'-e~ 2A F^ = 0, 
(e* Ftr)' = , 

h^ + (3A'h + h')$' a -m 2 J a = (3.7) 

with a constraint, 

* t fa - F tr ~ 2AV - 6h(A') 2 = -6+t *&& ■ (3-8) 

a=l a=l 

The third equation in (13.71 ) can be solved by 



* rt = 2qe- A{ - r \ (3.9) 



leading to the set of equations 



2 -i 



h" + 3A'h'-4q 2 e- 4A = 0, 

h^+(3A'h + h'W a -m 2 a Q> a = 0, (3.10) 

which will be the starting point of our subsequent analysis. 
The black brane solution, 

A(r) = r, h(r) = l-e- 3r+3r », F rt = ty a = 0, (3.11) 

describes the uncharged sector of the ABJM theory at finite temperature. Due to the conformal 
symmetry of the black brane background, the finite temperature phase here depends on only one 
dimensionful parameter which can be taken as the temperature T. Thus this uncharged sector 
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possesses only one finite temperature phase corresponding to the high temperature limit. The 
temperature T is identified with the Hawking temperature of the black brane, 

T=±h'(r H )e A ^ = ^-e r », (3.12) 
471 4% 

where r = r# is the location of horizon. The Bekenstein-Hawking entropy density is given 
by the expression S in (12.3b . The energy density, the free energy density and the pressure are 
related to the entropy density by £ = ^ST = 2p = —2[F as simply dictated by the conformal 
symmetry of the background. Some probes of system are investigated in Ref. j6l by studying 
the response of the scalar and current operators to the external perturbation. The static length 
scales including the true mass gap as well as the Debye mass are further studied in Ref. 0. 
Our main interest of this note is the RN black brane solution, 

A(r) = r, h(r) = 1 - £e~ 3r + q 2 e~ 4r , F rt = 2qe^\ §a = 0, (3.13) 

which is an exact solution of the original lOd equations of motion. The parameters £ and q here 
are respectively proportional to the mass and the charge densities of the RN black brane. The 
horizon is located at r = r# with h(rjj) = satisfying explicitly 



£e - 3rH = l+q 2 e- 4rH . (3.14) 



The minimum of h(r) occurs at 



e r,n = -r—!y , (3.15) 



3£ 

with h'(r m ) = 0. No nakedness condition for the mass and charge requires h(r m ) < leading 
to the inequality 

where the inequality is saturated at zero temperature. The Hawking temperature of the RN black 
brane becomes 

3ee~ 2rH / 4a 2 \ 



4% V 3£ 
The entropy, the energy and the charge densities read 



S = ^e 2r », £ = 4, p=% (3.18) 

K K K 

with J = £ — Ts and /jg = where /jq is proportional to the field theory chemical 

potential /j. Hence the no nakedness condition now takes a form 



K 2 £\ 2 /K 2 P\ 3 



In Ref. |fT8l , this character is attributed to build-up of a fermi surface for a finite number of 
the fermion number density. The weakly coupled massless fermions at finite temperature in 
general do satisfy an analogous inequality with precisely the same powers. Further support of 
the picture comes from the fact that the specific heat at low enough temperature is linear as 

Cv=JvT 1 (3.20) 

which is another important characteristic of the fermi surface. 

However, there are some additional properties which cannot be understood from the fermion 
picture. The entropy density even at zero temperature remains finite; namely, 

5(r = 0) = ^=p, (3.21) 

with a finite size of horizon at e 2r »( r=0 ) = -2=. Where this entropy comes from even at zero 
temperature is not clear. Recall further that the U(l) current is conserved only in sum of bosonic 
and fermionic contributions together. Thus at weak coupling the build-up fermi surface would 
not be possible if the bosons were in an unbroken phase and, hence, putting charges to the 
bosonic states were energetically favored. However the symmetric unbroken phase without 
condensation will be problematic as argued in the previous section for the weakly coupled case. 
Therefore the system should be at some unbroken phase with some bosonic condensate at least 
for the weakly coupled regime. In later sections, we would like to study the nature of phases 
occurring at the RN black holes, in the strongly coupled side, by the condensation of some 
operator expectation values. 

Near boundary regions of large r, general asymptotically AdS solutions should behave as 

A(r) = air + ao-\ , 

h(r) =h(oo) +h 3e - 3A -\ , 

F rt (r) = 2qe- A , (3.22) 

and the discussion of the boundary data for the scalar fields will be specified below. The entropy, 
energy, charge densities and the temperature have the expressions, 



K 2 ' K 2 /z(°o)' 

p = ?f, T=^)^=. (3.23) 
In this computation all the thermodynamic quantities are measured with respect to the boundary 



time \ fhj°°)t such that one may bring the boundary metric of (13.61 ) in the standard form. 
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4 Thermodynamic stability of the RN black brane 



In this section we shall discuss the thermodynamic stability of the RN black brane solution 
(13.131) . The thermodynamic stability is ensured if the Hessian (second-derivative) matrix of 
the energy with respect to its thermodynamic variables has no negative eigenvalues. With any 
negative eigenvalue, the system becomes thermodynamically unstable under small fluctuations 
that drive the system toward some other stable point lfT9l l20l |2T1 . In the type IIB theory, 50(6) 
R-charged black brane solution is available. But there it is observed that the R-charged black 
brane solution exhibits thermodynamic instabilities at a temperature lower than a certain critical 
value ifTTI . The fate of black brane in this unstable regime has not been known up to now. 

Unlike the case of this type IIB counterpart, we find that our RN black brane solution is 
thermodynamically stable. To show this, note first that the energy density can be expressed as 



3 



(271) 



2~2 



/ 7C z p 



3 V 1 c2 

2 v 



1 + (4.1) 



The components Hu of the Hessian matrix are given by 

9 2 £ 3k / % 2 p 2 

Hn = ^2 = „ , n , 3 i ( 1 + — — 



ds 2 4(2tu)M V 

2 

K7C 



a 2 £ 



(2tc)M' 



5 P 2 (2tc)M 



K 2 / K 2 p 2 



The determinant of Hjj then becomes 

together with tr H > 0. This proves that the two eigenvalues are positive definite. One may also 
consider turning on small quantities of the three further charges p a (a = 1 , 2, 3) in the Cartans of 
the SU (4) R-symmetry. Including these contributions to the quadratic order, the energy density 
has the expression, 

3 ? 

/ 7T 



,3 V 1 + 72(P +P«P«) • (4-4) 

Using now 5x5 Hessian matrix, one may easily verify that the system at p fl = is thermody- 
namically stable even in this enlarged space. Thus we conclude that the RN black brane solution 
is thermodynamically stable. 
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5 Geometrical instability of the RN black brane 



In this section, we shall investigate possible geometrical instabilities of the RN black brane sys- 
tem with a particular probe mode turned on. Depending on the mass of the supergravity mode, 
the RN black brane may become geometrically unstable below a certain critical temperature T c 
leading to a new black brane solution wearing nontrivial hairs. 

For scalar modes, arising of the instability may be understood as follows lfT5lfi6ll22~ll . Note 
that the usual geometrical stability condition for the AdS^+i spacetime is given by the so called 
Breitenlohner-Freedman (BF) bound [|23l . 

d 2 

-— < m 2 , (5.1) 

which is indeed respected by any scalar mode of the present theory. This is the relevant condi- 
tion for the stability of the near boundary region of the RN black brane solution that is asymp- 
totic to AdS4. On the other hand, the near horizon geometry of the extremal RN black brane in 
d + 1 dimensions is given by AdS2 x W l ~ l with a scaled AdS radius of 1/ y/d(d— 1). Hence 
for scalars, the BF bound of this region is violated if 

2 d(d-l) 

m < — — -. (5.2) 

4 

Therefore for our case of d = 3, the instability below a certain temperature may occur if the 
mass squared of a bulk scalar is ranged in 

9 i 3 
-- A <m 2 <--. (5.3) 

For the higher spin fields with spin s > I, one may show that there is no potential instability 
once they are neutral under the U{\) gauge field of the RN black brane. Thus scanning the 
supergravity modes in Tabled] one finds that the possible gravitational instabilities are limited 
to the case of scalar fields with m 2 = —2. As explained in Section 3, these bulk scalars are dual 
to the field theory operators of dimension A = 1 and 2. 

To show the instability of the black hole, let us study the positivity condition of the energy 
functional for the scalar field. The energy of a scalar field reads 



E = j drdxdyV^lW'li^+g'^^+g^Mf+g^XMf + 'nV 



(5.4) 



where dots and primes denote derivatives with respect to t and r respectively. If there exists any 
normalizable (probe) mode cp of the scalar field (|) which makes this energy functional negative, 
the geometrical instability of the RN background can be triggered driving the system to some 
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new stable configuration. In order to find a possible mode of instability, we take cp as a function 
of r only and look for the negative fluctuation mode satisfying 

(e 3A h(r) <p , y + 2e 3A <p-x§<p = 0, (5.5) 

where Ko is a real constant and we set m 2 = —2. Of course one may turn on the spatial fluctuation 
by considering the probe field depending on x and y by e l ( k * x+k y y > but this will only increase the 
energy of the system. Hence the above consideration will be sufficient. 

The boundary conditions are crucial for a determination of the solution of (15.51) . For r = °°, 
we note that the behavior of scalar fields in the near boundary region takes a form, 

cp ~ s A (x)e^ 3 - A ^ + OA (jc)e- AA ('') + • • • , (5.6) 

where ■ ■ ■ denote higher order terms. From the standard dictionary of AdS/CFT, the presence 
of sa corresponds to turning on an external source term for the dual operator Oa(x) while oa 
represents the operator expectation value (Oa(x)} Sa in the presence of the source sa(x). In our 
present problem, we would like to consider the system without introducing the source term 
and, thus, our choice of the boundary condition here and below is sa = at r = °°. The other 
boundary condition is for the horizon of black brane. Basically we shall require nonsingularity 
of the configuration there. For this, we need to expand the equation (15.51) into a generalized 
power series in r — r#, and then we can obtain boundary condition for the horizon as 

2 - }&e~ M ( rH } 

<™ = — 

<P"M = ° 7T 2 (3A\r H )h\r H )+h>\r H )+2-K 2 e- 3A ^)< ? (r H ) 

2{h'{r H )Y V ' 

\y \ ( 5 - 7 ) 

2h\r H ) 

by setting the coefficients of the negative powers of r — r# to zero. Note that the denominators 
are given by h\rn)'- If it is too small, we may have trouble in numerical computations. In order 
to avoid potential numerical instabilities, we use another form of RN solution given by 

A/ . ,„ 2 , l-(l+A-3A(r) + 2 g -4A(r) 
A r ) = (3-q 2 )r, h(r) = ^ ^- ^ . 

(3-r) 2 

2q 



A f = -3^( e - A W-l), Cb fl = 0. (5.8) 
This solution can be related to (13.131) by the coordinate transformation, 



r 



■ q 2 e- 4r n 



t> = e r »(3-q 2 e-^] 



x' = e rH x, y' = e rH y, (5.9) 
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with the redefinition of the parameters, 

q' = q e - 2rH , e' = Ee~ irH . (5.10) 

After the transformation, we drop primes for the notational simplicity. The temperature for this 
system reads 

T = ^(3-q 2 ), (5.11) 

and the location of the horizon is at r = 0. Here and below, we shall use this rescaled background 
for the numerical analysis. 

Now one can solve the probe equation (15.51) numerically to find the temperature below which 
the negative mode begins to develop. For the numerical analysis, we use the standard shooting 
method based on a Mathematica coding. The results are as follows: The on-set of instability 
occurs at 

r c = 0.0395(7) [? c = 1.582(0)] (5.12) 

for the A = 1 operator. For the A = 2 operators, we found 

T c = 0.0003(5) [q c = 1.7307(7)]. (5.13) 

The differences in the numbers of significant digits of T c and q c arise due to the fact that T c is 
proportional to 3 — q c . 



6 Phase structures and critical exponents 

In the previous section, we have established the instability of the RN black brane in the pres- 
ence of bulk scalars with m 2 = — 2. As temperature decreases below the critical temperature, 
the scalar fields begin to develop a nontrivial profile that affects the original RN black brane 
geometry. As we shall see later on, the boundary CFT undergoes a phase transition by the 
condensation of the expectation value of the dual field theory operator 0&(x). 

In this section, we shall investigate this phase transition by looking at the supergravity solu- 
tions. Since the changes in the thermodynamic quantities like the entropy, energy and so on are 
encoded in the geometry, the probe analysis of the previous section alone is not sufficient and 
inclusion of the full back-reaction to the geometry will be essential. Our study of the relevant 
solutions will be mainly based on numerical analysis. 

Our starting point of the analysis^ is the set of equations in (13.91) and (13.101) where we turn on 
only one scalar field (]) with m 2 = — 2. To set up the problem completely, one has to also specify 

3 We use here our original coordinates without performing the coordinate transformation in (15.91 and (15.10b . 
Then, when the scalar field vanishes, we shall directly obtain the exact RN solution ( 15.8b by fixing some freedoms 
in our coordinate choice. 
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the boundary conditions. Let us begin with the horizon side. At the horizon, h(r) has to be zero 
at some finite r = r# which is the coordinate location of the horizon. Using the translational 
freedom of the r coordinate, we shall choose 



r H = 0. (6.1) 

In addition, we note that one has the scaling freedom for x M = {t,x,y) and the coordinate r to 
generate a new solution. Using this we shall fix 

A(0) = and h'(0) = l. (6.2) 

At the horizon we basically require that all the fields in (13.101) should be regular. Setting (|)(0) = 
u, the regularity leads to 

2u 



f(0) 



h'(0) 



3_ a 2 -44(0) , u 2 



and 



h"(0) = -9-3u 2 + lq 2 e- AA ^\ 

A " (0) " W 

4>"(0) = -^(\+2q 2 e-^). (6.4) 

(Solving ( 13.91) and (13.101) with the coordinate and boundary conditions in (I6.1I) - (I6.4I) for the 
case of vanishing scalar field leads directly to the exact RN solution in (15.81) without need of the 
transformation in (15.91) and (15.101) .) At r = °°, we shall impose the behaviors of fields in (13.221) 
that are required for the asymptotically AdS spacetime. Again note that the scalar field in the 
near boundary region takes the form 

4> ~ * A (*)<T (3 - A)A(r) + 0A(*)*" M(r) + " " • , (6-5) 

where • ■ ■ denote higher order terms. (The interpretation of s&(x) and o^{x) in the boundary 
CFT is the same as that in the previous section.) We shall set the source term s& (x) = 0, which 
corresponds to the boundary system without an external source term. By this last condition, u 
will be determined as a function of q by u = u(q). 

Now adopting the shooting method based on a Mathematica coding, we perform a numerical 
analysis for A = 1 and 2 cases separately. For each case, we find one set of solutions that 
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is parameterized by the value of q. The resulting functions u(q) in (u, q) plane are drawn in 
Fig- HI F° r eacn case, there is a critical value of q c beyond which the corresponding scalar field 
begins to develop a nontrivial profile. Namely if q < q c , the exact RN solution in (15.81) . which 
satisfies all the coordinate and boundary conditions in (I6.1I) - (I6.4I) . remains intact. This part is 
depicted in Fig.l by the vertical blue solid line for q < q c . If q > q c , the RN black brane is 
modified by wearing a nontrivial scalar hair with u(q) ^ 0; In Fig.l, the blue dots represent the 
data set which we obtained by the numerical analysis. The blue solid curve represents the fitting 
function that is obtained by the standard curve fitting method in Mathematica. 




Figure 1: The functions u{q) in (u, q) plane are depicted in the left and the right sides respectively 
for the A = 1 and the A = 2 cases. For each case, the development of the scalar profile is represented 
by nonvanishing u{q) for q > q c ; The dots represent our numerical data set while the solid curve is for 
the fitting function obtained by the standard curve fitting method in Mathematica. Below q c , u{q) = is 
indicated by the blue solid line, which corresponds to the RN black holes in (15.81) . In addition, we marked 
the numerical values of q c by the red circles, which are given by 1.581(6) and 1.7308(3) respectively for 
the A = 1 and the A = 2 cases. 



Thus the critical temperature may be evaluated by using the RN black brane solution (15.81) 
with q = q c leading to 

T c = ^(3-ql). (6.6) 
By the standard curve fitting method, we found 

T c = 0.0396(4) [q c = 1.581(6)], (6.7) 
for the A = 1 scalar field, and 

T c = 0.0003(4) [? c = 1.7308(3)] (6.8) 
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for the A = 2 scalar fields. In principle, these critical values should agree precisely with those 
from the probe analysis of the previous section. Therefore the numerical precision of our analy- 
sis can be estimated by comparison of the numerical values from the two methods. We see that 
the A = 2 critical temperature, whose value is closer to zero, has a poorer numerical accuracy. 

In order to show further numerics, instead of using the charge density p = 2q/K 2 , we shall 
use simply q which differs from the actual charge density by a constant factor 2/k 2 . (See (13.51 ) 
for the definition of K.) In the above sets of solutions, both q (or the charge density) and the 
temperature T = T(q) changes as one changes q along the transition. But what we want is to 
fix q (or the charge density) while changing the temperature of the system along the transition. 
In order to generate such sets of solutions, we use the coordinate transformation of the gravity 
system by 

-> (6.9) 
Note that, by the scale transformation, the thermodynamic quantities transform as 

T -> l s T, q -> ijq, £ -)■ X^E, 

S -> X 2 5, o A -> A^oa, sa -> kl~ A SA, (6.10) 



and, hence, 



(6.11) 

Ota dSA 



Using this scaling transformation, we shall fix q = 1 by choosing X s = 1/ ^fq. ( The quantities 
carrying a tilde are the ones after the scale transformation (16.101) . ) The new sets of scaled 
solutions are now parameterized by the rescaled temperature 

f = 4= (6.12) 

with fixed charge density p = 2/k 2 . Using the values in (16.71) . we found the rescaled critical 
temperatures f c = T c / ^fq~ c as 

T c = 0.0315(6) (6.13) 
for the A = 1 scalar field. For the A = 2 scalar fields, we found 

T c = 0.0002(5) (6.14) 

using the values in (16.81) . 
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Figure 2: The phase diagrams in the left and the right sides are respectively for the A = 1 and the A = 2 
scalars. The system undergoes a symmetry breaking transition from a symmetric phase to a broken phase 
as the temperature is lowered below the critical temperature. For each case, the blue dots represent our 
numerical data set and the blue solid curve depicts the fitting function that is obtained by the standard 
curve fitting method in Mathematica. We marked the numerical values of t c by the red circles, which are 
given by 0.315(6) and 0.0002(5) [see (16.131 ) and (16.141 )1 respectively for the A = 1 and the A = 2 cases. 



Below the critical temperature, one finds a development of expectation value oa which sig- 
nals a phase transition. As we shall see details later on, this basically corresponds to the spon- 
taneous symmetry breaking transition of the SU (4) R-symmetry in which the condensate oa 
plays the role of the order parameter. In terms of the scaled variable oa = o^/q^l 2 , the phase 
diagrams are depicted in Fig. |2] for the A = 1 and the A = 2 cases. 

As expected, the transition for A = 1 occurs at a higher critical temperature. Then the 
transition for A = 2 cannot be treated separately. Namely one has to turn on both of the A = 
1 , 2 scalar fields around the region of the A = 2 transition in which the A = 1 scalar field has 
developed some finite amount of profile representing the condensation. However remember 
that the scalar part of our starting Lagrangian is only valid up to quadratic orders. Hence our 
treatment loses its validity around the region of the A = 2 transition. Though it is an interesting 
problem to clarify further, we shall leave it to the future investigation. For the remaining we 
shall be focusing on the transition involving the A = 1 condensation. 

The phase transition is of second order as in the usual cases of symmetry breaking transition. 
We note that natures of phase transitions are in general classified by their critical exponents. 
Here we compute numerically the exponents a, P and y respectively defined by 



\T-Z 



dT 

o\ ~ \T — T C 



c 



do i 
Xl = 3— 

OSl 



.51=0 



\t-t c \-i. (6.15) 
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For the numerical estimation of the exponents a, (3 and y, we use the numerical data sets 
respectively of the forms (log(T — T c ), \og(dE/dT)), (log(T — T c ), logoi) and (log(T — T c ), 
logxi). Using the fitting with linear least squares, we estimated the relevant slopes as well as 
their standard deviations. From the corresponding data set for o\ in Fig. |2], one finds 

% =0.04978 ±0.0032, (6.16) 

which is in a good agreement with the mean field value P = 1/2. Here and below the subscript n 
indicates that the relevant exponent is obtained by the numerical analysis. The left side of Fig. |3] 
shows the behavior of the energy as a function of temperature in the vicinity of the transition 
region. Form the corresponding data set, the exponent for the specific heat is identified as 

a n = 0.012±0.018, (6.17) 

which agrees well with the mean field value a = 0. The right hand side of Fig. [3] shows the 
temperature dependence of the susceptibility log %i with respect to the logarithm of temperature. 
Again from the corresponding numerical data set, one finds 

y n = 1.011 ±0.015, (6.18) 

which agrees well with the mean field value y = 1 . One can check that 

(a±2p±y), 3 = 2.019 ±0.024 (6.19) 

which is consistent with the so called Rushbrooke scaling law 

a±2(3±y=2. (6.20) 

Finally we shall reconfirm the above exponents by the analytic treatment. We note first that, 
due to the boundary condition (|)(0) = u, the scalar field behaves as <^> <=< u when u < 1, Then 
inspecting the equations of motion (13.101) together with the boundary conditions (16.31) and (16.41) . 
one finds that si and o\ should behave as 

s\ = u (oo(q) + a 2 (q)u 2 + 0(u 3 )) 

o x = u(b (q)+b 2 (q)u 2 + O(u 3 )), (6.21) 

before imposing the last boundary condition of s\ =0. For \u\ <C 1, the coefficients of the u 3 
terms in si and o\ are basically controlled by the u 2 terms of A'(0) and h'(0) in (16.31) and (16.21) . 
which one may argue by a careful examination of the equations of motion in (13.101) : By this 
consideration, one may show that 

a 2 (q)=a 2 o + 0(q-q c ), b 2 (q) = b 2 o + 0(q-q c ) (6.22) 
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Figure 3: The left figure shows the scaled energy density as a function of temperature. The right hand 
side shows log \ *Jq%\ \ in ( 16.151 ) with respect to log (j / yfq — T c /^/q^j. For each case, the blue dots 
represent our numerical data set and the blue solid curve depicts the fitting function that is obtained by 
the standard curve fitting method in Mathematica. 

with <32o ^ and &20 0- Then, since s\ = has a nontrivial solution u(q) ^ only for q > q c , 
ao(q) should behave as 

a (q) = a i (q - q c ) + 0[ (q - q c ) 2 ] (6.23) 

with floi a 20 < for \q — q c \ 1 . (Here we assume the existence of q c , which is justified by our 
numerical analysis.) Since we know that the o\ = boundary condition instead of si =0 leads 
to a different value for the critical charge, bo(q) has the expansion of the form 



£o(<?) =boo + 0(q-q c ) 
with boo 7^ 0- This argument shows that the si = condition leads to a solution 



-QQl 
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(q-q c )* 0- + O(q-q c )) q > q c 
q < q c . 



Thus, 



\T-T C \^ 



(6.24) 



(6.25) 



(6.26) 



which implies that P = i. Since the scalar field contribution to the energy density is of order 
u 2 , the energy density has to be of the form 



£ = eo + e\ (q - q c ) + e2U 2 + ■ 



(6.27) 



where the first two terms are from the original RN black brane with u = 0. We then conclude that 
a = since the leading power of specific heat is zero. By changing u = u(q) to u = u(q) + du, 
si and o\ in (16.211) vary by 

8^i = 3a2o(u(q)) 2 du 1 801=^008" (6.28) 
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to the leading orders. Then the susceptibility behaves as 

80 1 _ frpp 1 ^ ir-T^r 1 (6 29) 

851 3a 2 o (u(q)) 2 

which implies that y = 1 . This proves that our exponents are those of the mean field theory. 



7 Interpretations and discussions 

In the previous section, we have established the phase transition whose exponents belong to the 
universality class of the mean field theory. For the resulting expectation value of A = 1 operator, 
let us introduce a notation 

Mu{x) = (O'Ax)) (7.1) 

where Oj is given in (13.21) . Since Mu = Mfj and Mu = 0, the 4x4 matrix M is Hermitian and 
traceless. We have seen that the boundary CFT undergoes a phase transition. Above the critical 
temperature T c , the condensate is vanishing, i.e. M = while at lower temperature, M ^ 0. 
This is the spontaneous symmetry breaking phase transition where the SU (4) R-symmetry of 
the ABJM theory is broken by the presence of the condensate M . 

Using M as an order parameter, the phase transition may be effectively described by the 
Landau free energy F 

F/T = f d 2 x[ TrVfM" • VfW -jlxMrt + 2r (T -T c )TrM 2 + go(TrM 2 ) 2 ^) , (7.2) 

where x = (x, y) and the 4x4 traceless Hermitian matrix is for the external source term in the 
adjoint representation of SU (4) . The form of the quartic term is determined by the symmetry. 
Due to the symmetry of the gravity solution, it should be invariant under the transformation of 
M by the SU(4) R-symmetry. The Landau free energy is minimized if 



T > T c 

(7.3) 

^{T c -T)2fi T < T c 



where h is a 4 x 4 traceless Hermitian matrix with Tr h 2 = 1 . This implies that P = 1 . Noting 



Cy = -rg| With 



[T-T c ) 2 , (7.4) 



one has a = 0. Similarly one can check that y = 1 from the definition of the susceptibility. One 
can further compute the correlation length scale ^, 



£=(2r (r c -r))-2 ~ |r-r c |" v (for T < T c ) , (7.5) 
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from which one finds v = 1/2. It will be interesting to test the prediction of the mean field theory 
by the direct analysis of the bulk gravity system. But on top of the new background obtained 
from the numerical analysis, further numerical analysis of solving the scalar fluctuation equation 
has to be performed. However this requires some improvement of the current method to reach 
the required numerical precision. 

As we described in Section 3, our U(l) charged black brane describes a system where the 
t/(l) number is carried by both bosonic and fermionic degrees. The build-up of the fermi 
surface at weak coupling side would not be possible if the bosonic degrees were in a unbroken 
phase. Indeed we have argued that there will be a condensation of the elementary degrees at 
weak coupling with low enough temperature. For the strong coupling side, we know that the 
basic degrees that are weakly coupled among themselves are now organized by the supergravity 
modes. We found that the symmetric phase of the RN black brane again becomes unstable 
when T < T c . The effective mass squared of the Landau free energy description becomes 
negative in this region driving the symmetry breaking phase transition by the condensation of 
the operator expectation value. This partially explains the fate of the bosonic U( 1) charged state 
of the ABJM theory at strong coupling. But as observed before, there is a finite entropy of the 
extremal RN black brane at zero temperature. However due to the phase transition, ordering of 
degrees will lead to a reduction of entropy in general. Hence there can be a chance of having a 
zero-entropy system at zero temperature. One interesting point of study is then the modification 
of the RN black brane by the scalar condensation in the near zero temperature region. Since 
the scalar part of our starting Lagrangian is only valid up to quadratic orders, we cannot answer 
this question at the moment. One has to improve our understanding of higher-order scalar 
contribution to the Lagrangian first. 

Another interesting aspect is the build-up of the fermi surface as mentioned before. The 
picture presented in Ref. [Tl8l is seemingly plausible, which mostly concerns the region of near 
zero temperature, and works only for a pure RN black brane system. But now close to zero 
temperature, our black brane solution is significantly modified by the condensation of nontrivial 
scalar profile. Further studies of the low temperature region appear interesting, in particular with 
focus on the question whether the Fermi picture is still valid. 

Finally we comment on the theorem by Coleman-Mermin- Wagner- Hohenberg E4l . The 
theorem states that a continuous global symmetry cannot be broken spontaneously in 1+1 di- 
mensions at zero temperature and 2+1 dimensions at finite temperature. Our example appears 
to be contradicting with what the above theorem states. However, we are working in the strict 
large N limit, which makes the dual gravity system completely classical. This violates some 
assumptions of the above theorem as pointed out in Ref. Il25l . Our example here is reminiscent 
of the clash of the unitarity in the explicit time dependent black hole solution which describes a 
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thermalization of a boundary field theory ||26l . There again the trouble stems from the large N 
limit of the boundary field theory. 
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